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The  linear  properties  and  saturation  mechanisms  of  the  lower-hybrid  drift 
instability  have  been  investigated  using  a  one-dimensional  particle-hybrid  simu¬ 
lation.  For  low  drift  velocities  (  \\,«v„  ),  ion  trapping  and  current  relaxation 
(  i',,— 0>  are  competing  processes  for  stabilization.  If  the  relative  electron-ion 
drift  velocity  is  kept  constant  in  time,  ion  trapping  causes  the  end  of  instability 
growth;  if  this  drift  is  allowed  to  decrease  consistent  with  momentum  balance, 
then  saturation  is  due  to  current  relaxation. 
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LOWER-HYBRID  DRIFT  INSTABILITY  SATURATION 


MECHANISMS  IN  ONE-DIMENSIONAL  SIMULATIONS 

Yu-Juian  Chen  and  C.K.  Birdsall 

Electronics  Research  Laboratory 
University  of  California,  Berkeley.  C A.  94720 


I.  INTRODUCTION 

The  lower-hybrid  drift  instability1  may  play  an  important  role  in  the  dynamics  of  plasma 
experiments,  such  as  fl-pinches  and  field-reversed  configurations.  According  to  the  linear 
theory  and  experiments,  during  the  implosion  phase  of  a  theta-pinch  this  instability  is  responsi¬ 
ble  for  the  (lute  observed  in  the  current  sheath  region’"4  .  Two  dimensional  particle  simula¬ 
tions  of  this  instability  with  a  finite  plasma  beta  value  for  large  drift  velocities  were  performed 
by  Winske  and  Liewer'  .  During  the  post-implosion  phase,  it  is  considered  to  be  the  dominant 
instability"  .  If  this  instability  saturates  at  high  levels,  it  will  contribute  to  anomalous  tran¬ 
sport.  and  influence  the  resistivity,  for  example.  Saturation  mechanisms  of  the  lower-hybrid 
drift  instability  during  the  post-implosion  phase  is  the  main  concern  of  this  paper. 

In  the  post-implosion  region,  the  difference  of  the  cross-field  velocity,  vd  =  vdl.  —  vd,  .  is 
less  than  the  ion  thermal  velocity,  i\  =  ( T/ni I1  ;  .  Anomalous  transport  coefficients  in  this 
region  were  derived  from  quasilincar  theory”1  .  For  low-beta  plasmas,  Davidson"  predicted  that 
current  relaxation  (  v,  —  0.  )  and  plateau  formation  (  6/,/dv,  —  0.  )  in  the  ion  velocity  distri¬ 
bution  are  generally  competing  processes  for  stabilization.  If  the  initial  drift  velocity  >’,/  is  very 
small,  then  it  is  energetically  favorable  for  stabilization  to  occur  through  formation  of  the  pla¬ 
teau.  Otherwise,  plateau  formation  will  not  have  been  completed  when  the  cross-field  current 
relaxes  to  zero.  Electron  resonance  broadening  has  also  been  considered  as  a  possible  satura- 
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lion  mechanism  for  high-beta  plasmas'* 1,1 . 

In  this  paper,  one-dimensional  particle-hybrid  simulation  studies  of  the  lower-hybrid  drill 
instability  are  presented.  Our  simulation  model  is  a  slab  with  a  constant  density  gradient;  the 
ions  are  unmagnetized  particles,  shielded  by  strongly  magnetized  electrons  through  the  linear 
electron  charge  density  susceptibility,  x, .  Ions  are  initially  in  a  steady  equilibrium  stale  with 
the  ion  diamagnetic  drift  velocity  cancelled  by  the  L*B  drift,  corresponding  to  electrostatically 
confined  ions.  At  small  amplitudes,  our  simulations  show  good  agreement  with  linear  theory, 
such  as  the  linear  growth  rate,  the  mode  frequency,  and  the  influence  of  finite  beta  effects  asso¬ 
ciated  with  the  nonresonant  electron  orbit  modifications"  15  .  Ignoring  effects  of  finite 
beta,  it  is  found  that  the  end  of  wave  growth  is  due  to  ion,  trapping  when  both  the  density  gra¬ 
dient  and  the  cross-field  current  velocity  r,,  are  kept  constant  in  time.  When  \\,  is  allowed  to 
vary  in  time,  our  simulations  show  that  stabilization  of  the  wave  occurs  through  current  relaxa¬ 
tion. 

In  Section  II.  A,  a  brief  review  of  the  linear  theory  of  the  lower-hybrid  drift  instability  is 
given.  A  description  of  the  simulation  method  is  presented  in  Section  II. B.  We  solve  the  linear 
dispersion  relation  of  our  simulation  model  (  with  a  multibeaming  system  )  numerically.  Com¬ 
parisons  of  observed  linear  properties  of  the  growing  waves  with  linear  theory  are  made  and 
given  in  Section  111.  Investigation  of  various  saturation  processes  is  presented  in  Section  IV. 
At  low  drift  velocities,  ion  trapping  causes  the  end  of  wave  growth  when  drift  velocities  are 
kept  constant  in  time.  It  is  found  that  the  saturation  level  predicted  by  quasilinear  theory  is  the 
saturation  level  of  the  most  unstable  mode  for  ion  trapping.  A  detailed  explanation  is  given  in 
Section  III. A.  An  analytic  formula  for  the  saturated  field  energy  spectrum  in  k  space  for  ion 
trapping  when  »•,/«»•,  is  derived.  Section  III.B  is  devoted  to  a  description  of  how  current 
relaxation  can  stabilize  the  instability.  The  simulation  method  for  allowing  r,  to  vary  in  lime  is 
also  described  Simulated  saturation  levels  agree  fairly  well  with  the  predictions  of  quasilinear 
theory  .  Conclusions  are  given  in  Section  V. 


II.  I.INKAR  TUTOR*  AND  SI  Ml  I  A  I  ION  MOI)KI. 


A.  Review  of  linear  Theory 

Our  model  is  a  slab  with  a  constant  density  gradient  in  the  y  direction  (I  ig.  I )  The  elec¬ 
tron  drift  across  the  magnetic  field  B„  :  with  mean  fluid  velocity 

*'/,  =  >7  +  i-,  <11 

where  17  is  the  l:*B  drift  velocity,  and 

v.  = - - —  ~ I n ( //,  )  (2) 

in,  w  ,  or 

is  the  electron  diamagnetic  drift  velocity,  w,,  is  the  electron  cyclotron  frequency,  in,  is  the  elec¬ 
tron  mass.  /  is  the  electron  temperature  and  n (y )  =  n,~n  is  the  density.  There  is  no  net 
ion  current  in  the  \  direction,  i.e., 

i\  =  17  +  i  -  =().,  13) 

which  is  a  consequence  of  the  ion  pressure  force  being  balanced  by  the  ambipolar  electric  force 

in  the  y  direction,  r-  is  the  ion  diamagnetic  drift  velocity  as  given 

__L.A,n (4) 

m  w.  dr 

where  l  .  in  .  w  are  the  ion  temperature,  mass  and  cyclotron  frequency,  respectively.  The 

electrostatic  simulations  were  performed  with  E\  and  displacement  only  in  the  x  direction. 

The  lower-hybrid  drift  instability  is  characterized  by  strongly  magnetized  cold  electrons 

and  unmagneti/ed  hot  ions  with 

I  I  , 

(fi  «  j  o)  4-  r y  \  <<.  to  , 

A  'll  '  »  I  ,  k'tl.  ’  S*  I.  <5» 

/  »  /;. 

With  these  inequalities,  the  dielectric  function  for  cold  electrons  and  Maxwellian  ions  is 
expressed  as 

/)  t  A  ,<o)  ~  l)f,'  ( A  .to)  -t-  //)/ 1 A  ,o>) 


(6) 


The  real  pari  of  the  frequency  is  determined  to  zeroth  order  in  «  1  by 


<o,„  I  to, 

1  +  -4-  +  - —  =0 

u>  k  ‘A  j)l  (u  ,—k  \  f 


The  solution  is 


to,  '  ^  rAl’/  —  to  , 

A:+A,;, 

The  growth  rate,  y=-Di/(Bl)H/di o)„  ,  is  given  by 


-V? 


A  /  A'm  ^  .  l'/ 

—  .  i.i  '  "7  ' *  >  a,. 


2  (1-U-7a,;>’  I  A,„  v, 


,  ■>  1  1  :  /  > 

A,r,  -  — t--  ,  ,  2  <*»«./  *'/. 

A/>,  I +to/,,./<o,,. 

is  the  wave  number  of  the  most  unstable  mode,  and 


is  the  lower  hybrid  frequency. 


B.  Simulation  Method  and  Linear  Dispersion 

In  our  simulation,  the  ions  are  unmagnetized,  fully  nonlinear,  and  modeled  by  particles. 
The  relative  electron-ion  drift  velocity  is  >■,,  =  r(/,  —  r(/(  =  r,  when  T,  »  T,  .  Electrons 
were  represented  as  a  warm  fluid  with  linear  susceptibility  (A, to)  which  is  given  by 


<A,to)  = 


I  -  IJb)e~h  to,;.  A-i-v(A) 


A -2v;  to  —  A- »7 


from  the  local  theory,  where 


vy(fc)  = - ‘—IJb)v-h  <141 

/«,  to, ,. 

a  a  /,<*)  a  /,(/») , 

x  f-ln«  -  l-Ml-y— r)  -  f-|nr,  />(i- t-—)  . 

a.c  a.r  /„(M  ay  /„(/>) 

I, Ah)  is  the  modified  Bessel  function  of  the  first  kind  of  order  n.  The  Id  Fourier  transformed 


Poisson  equation. 
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1  +  X,  <A.to)  </>*(/)  =  ~~rf*k 


may  be  written  as 
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by  replacing 

to  with 

t  9 

the  operator 

ai 

kv/ 


d>k  (t) 


(15) 


(16) 


V  . 

%  -  k'l: 


dimensional  hybrid  simulations  for  the  drift-  cyclotron  instability  have  been  carried  out  by 
Cohen  et  al.IJ  I' . 


Quasilinear  theory7  *  predicts  that  the  field  fluctuations  of  the  lower-hybrid  drift  instability 
saturate  at  a  very  low  level  compared  with  the  thermal  energy  for  r(/ « i’„  .  In  order  to 
observe  growth  in  the  linear  regime  in  our  simulation,  it  was  decided  to  load  particles  orderly  in 
phase  space  (quiet  start),  with  a  Maxwellian  distribution,  in  order  to  reduce  the  initial  fluctua¬ 
tion  level.  However,  the  quiet  start  Maxwellian  consists  of  many  discrete  beams  which  can  pro¬ 
duce  multibeam  instability  as  pointed  out  by  Dawson lf>:  hence,  we  needed  to  analyze  the  nor¬ 
mal  modes,  of  the  ion  beams  and  the  cold  electrons  with  x<  from  Eq.  (13). 

The  ion  susceptibility,  x,  ,  is  obtained  from  applying  the  linearized  fluid  equation  of 
motion  and  equation  of  continuity  to  each  ion  beam  of  zero  order  drift  velocity  V,  and  density 
V  ,  with  the  well  known  result,  for  one  ion  beam 


4 

ni,  (to— A  ) , )’ 


(17) 


The  dielectric  function  is  DlA.to)  =  1  +  x,  +  x,  so  that  the  dispersion  relation  D  =  0  for  many 


ion  beam  is. 


l+X, 


v  477  A>~  1 

",  m,  (a,  ~  k  V , ) : 


118) 


The  dispersion  relation  was  solved  numerically  by  using  the  solver  supplied  by  Au-Yeung 
and  Friedman1  .  It  is  found  that  there  are  two  kinds  of  modes  in  the  system:  the  lower-hybrid 
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drill  wave  with  dispersion  affected  very  little  by  the  discrete  beams;  and  the  multibeam  instabil¬ 
ity  with  growth  rate  comparable  to  that  of  the  lower-hybrid  drill  instability.  In  Fig.  2  dispersion 
curves  are  plotted  for  the  lower-  hybrid  drill  mode  and  the  fastest  growing  mullibeaming  mode 
with  parameters: 

X  -  Xh  =  16384  . 

in /in,  —  1600.. 

io'Jio;,  =  I, 

A/i./A.v  =1.44  , 

»■/  / 1;,  =  (1.424263. 

r  =  o . 

I.„i  L/i  =  O.and 
=0  . 

L„  =  l/(91nii/9,i  ),  Ln  =  i/(9in B„/dy)  ,  and  L,  =  l/OInT/dy)  are  the  scale  lengths  of  the 
density,  magnetic  field  and  temperature,  respectively.  Sh  is  the  number  of  beams,  which  was 
taken  as  the  number  of  particles  because  in  the  quiet  start  Maxwellian  loading,  each  particle  was 
given  a  different  velocity,  in  order  to  minimize  the  multi-beaming  instability. 

III.  SIMULATION  RESULTS  AT  LINEAR  STAGE 

The  corresponding  simulation  results  for  the  mode  number  M  =  5  are  shown  in  Fig.  3 
(the  mode  energy  history  plot)  and  Fig.  4  (the  power  spectrum  vs.  frequency  plot  at  various 
times),  where  the  mode  number  M  is  defined  as  M  =  kL/  2n  .  Typical  parameters  used  are  the 
number  of  spatial  grids  N(i  =  64,  the  lime  step  w,„A  /'  =  0.2  ,  the  length  L  =  44.43  A„,  .  All 
the  results  presented  in  this  section  are  from  single  mode  simulation,  i.e.,  only  one  Fourier 
component  of  electric  potential  was  used  to  push  ion  particles.  Figure  4  shows  a  well  defined 
lower-hybrid  drift  frequency  at  w/a>(„  =  0.136  as  well  as  many  multibeam  mode  frequencies. 
As  time  goes  on.  the  lower-hybrid  drift  wave  grows  faster  than  the  mullibeam  modes  do,  and 
eventually  becomes  dominant.  We  measured  the  lower-hybrid  drift  mode  peak  values  at 
different  times,  and  calculated  its  linear  growth  rale.  The  simulation  frequencies  and  growth 
rates  are  shown  in  Fig.  2.  they  agree  with  the  prediction  very  well. 
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In  general,  the  growth  rate  for  a  given  mode  number  M  increases  as  17  / »>,,  increases.  Fig¬ 
ure  5  presents  mode  energy  history  plots  for  (a)  17/1,,  =0.57  and  (bl  0.85,  respectively.  The 
growth  rates  can  he  easily  measured  from  the  slopes  of  Figs.  5a  and  5b.  The  dependence  of  the 
linear  mode  frequency  observed  in  the  simulation  on  17  / 1,  and  comparison  of  these  results 
with  linear  theory  (Fq.  18)  for  the  mode  number  M  =  5  are  given  in  Fig.  6.  Mode  5  is  approxi¬ 
mately  the  most  unstable  mode  lor  the  parameters  we  used. 

The  addition  of  Unite  beta  alters  the  growth  rates.  The  beta  value  is  given  by 
/3  =  -21.  II. h  as  /  =  0.  In  Fig.  7.  dispersion  curves  are  plotted  for  17/17  =1.1314.  <<a) 
/.  /./(  =  0  .  and  (b)  /.„//.«  =  -0.40  where  the  negative  sign  means  the  gradients  of  density 
and  magnetic  field  have  opposite  directions.  Figure  8  gives  dispersion  curves  for 
i  :  .•  v  =  1.5556  .  (a*  /.,.//.«=  0  .  and  (b)  /.„//_#  =  —0.36.  Figs.  7  and  8  show  that  the  lower- 
hybrid  drift  waves  at  small  mode  number  M  are  stabilized  by  the  finite  beta  effects  associated 
with  the  nonresonant  V#„  electron  orbit  modification  (the  term-proportional  to 
(1  H  ){8B.iby)  in  Fq.  (6)  )  as  predicted  hv  Davidson  et  al."  “  .  Note  that  except  in  the  limit 
of  /  /  /  —0.  ,  they  predicted  that  there  is  a  critical  value  fi.-  for  the  local  plasma  beta  such  that 
the  lower-hybrid  drift  instability  is  completely  stabilized  for  fi>fi 

l\.  SAIT  RATION  MECHANISMS 

We  have  studied  several  saturation  mechanisms  of  the  lower-hybrid  drift  instability  in  an 
uniform  magnetic  field.  Simulation  results  show  that  the  instability  is  stabilized  by  ion  trapping 
when  both  the  density  gradient  and  the  relative  electron-ion  drift  velocity  r,  are  kept  constant 
m  time  because  of  external  circuits  for  example.  It  is  found  that  the  saturation  level  predicted 
by  quasilinear  theory  *  is  the  saturation  level  of  the  most  unstable  mode  via  ion  trapping. 
Allowing  the  electron  cross-field  drift  velocity  i-,/(  to  vary  in  time,  the  wave  stops  growing  at  a 


much  lower  level  as  the  current  relaxes 


A.  Ion  Trapping 


In  our  simulations,  the  potential  </<k(t  I  was  advanced  by  solving  Fq.  1)6).  Constant 
values  of  i,  ,  which  drive  the  lower-hybrid  drift  instability,  were  studied  first  In  single  mode 
simulation  for  the  most  unstable  modes,  there  are  oscillations  of  the  wave  energy  at  the  trap¬ 
ping  frequence  to  tl  ig.  5a  and  5b)  just  alter  the  time  saturation  and  a  vortex-like  structure 
emerges  in  the  ion  phase  space  plots  after  saturation  (Fig.  9a  and  9b)  for  1 7  /  v  =0.57  and  0.85. 
The  end  of  wave  growth  was  accomplished  by  ion  trapping.  All  of  the  phase  space  pictures  are 
presented  in  the  wave  frame,  v  —  >■ i  —  i<ni\iaiii .  where  r,,„  was  calculated  from  linear  theory.  It 
is  noted  that  vortices  in  Fig.  9a  and  9b  at  small  velocities  are  not  due  to  the  multibeam  instabil¬ 
ity.  because  those  vortices  would  appear  at  large  velocities. 

We  compared  these  saturation  levels  with  the  quasilinear  theory  as  given  by  Davidson" 
under  the  assumption  that  the  spectrum  is  strongly  peaked  about  the  fastest  growing  mode,  for 
v,  <  r.  Fig.  10a  shows  that  the  simulation  data  are  in  agreement  with  the  quasilinear  saturation 
level  e,  given  by  his  Fq.  (48) 


45v'jt 


s/I  I 


nl 

I  +<U/,  ij)~ 


(19) 


There  is  a  factor  of  v2  difference  in  r  in  our  version  of  Davidson's  equations  because  of  our 
differing  definition  of  T .  /  —m  iy 

It  is  necessary  to  explain  why  the  simulations  shows  that  nonlinear  saturation  was  due  to 
ion  trapping  but  the  saturation  levels  agreed  with  Davidsons  quasilinear  theory.  In  deriving 
the  saturation  level  <  Fq.  19).  Davidson  began  with  an  energy  conservation  equation  and  the 
only  real  invocation  of  quasilinear  theory  seems  to  be  the  specification  that  saturation  occurs 
when  the  distribution  function  has  been  "flattened"  around  the  mode  phase  velocity.  Such 
flattening  could  in  principle  be  due  to  a  variety  of  causes  besides  the  usual  quasilinear  diffusion, 
for  example,  trapping,  liquating  the  field  energy  to  the  ion  kinetic  energy  change  due  to 
"flattening"  gives  his  Fq.  (46). 


(I- 


J~  )  f* (/A  ( I  +  -  )e  ■  ( /— 00 )  =  - =r 

or  J  A,:.  45  s/ jt 


V  ?  1 


III 


(201 
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lo  calculate  i he  lot!  Mile  ot  the  above  equation,  |V>-  J.^n  assumed  that  the  spectrum  is 
sullicientlv  peaked  ;ihou  1  A  ■  = '  ';.ie  wavenumber  corresponding  to  maximum  growth  tor  the 

mm  1  vi|uilibiium  conditions)  that  l+A7A„-.  =  2  is  a  good  approximation  in  the  integrand  in  Fq. 
(20i.  Ibis  was  equivalent  to  a  single  mode  assumption.  Kq.  (19)  was  then  obtained 

Ihe  trapping  frequency  is  given  by 


iu 


w  here 


(21 ) 


,/,  = 


87T€  . 

A,;  1 


is  the  electric  potential  l 'sing  Fq.  (19),  we  obtain 


4 

1  A 

•7 

45v  77 

y/2V;i 

Using  y  from  Fq.  (9)  at  A=A.„.  Hq.  (22)  can  be  rewritten  as 


W 

_  4  1 

1  1 

_  8  _ 

M^V, 

y  ■ 

45\-  77  | 

v'2rr 

>7 

(22) 


(22) 


i 


1.975 


l  *■/  I 

When  i  /  <  r  .  the  trapping  frequency  is  larger  than  the  growth  rale  and  hence  larger  (hen  the 
bandwidth  Ju«=--y  as  well  Therefore,  ion  trapping  will  be  the  saturation  mechanism  when  the 
fastest  growing  mode  is  dominant,  and  e,  in  Kq  (19)  will  be  (he  saturation  level  due  to  ion 
trapping.  Comparison  of  l;q.  (22)  with  the  observed  trapping  frequencies  (as  in  Fig.  Ill  indi¬ 
cates  fairly  good  agreement  for  cases  in  which  17  <  t.  . 


Finally,  multi-mode  simulations  have  been  performed,  i.e..  all  possible  modes  are  excited 
at  t=().  'typically,  parameters  used  were  the  same  as  those  for  the  single-mode  run  The 
fastest  growing  mode  tended  to  reach  the  same  saturation  level  no  matter  whether  only  a  single 
mode  was  kept  or  all  modes  were  included  in  the  simulations.  In  the  wave  frame  of  the  most 
unstable  mode,  the  ion  phase  space  plot  has  a  vortex  formation  about  1  \  =0  after  saturation  for 
1//1  -0.57  as  shown  in  I  ig.  12a.  I  he  length  of  the  plasma  was  l.  =  2n  .  and  the  number  ot 


■  ]() . 


grids  was  Nti=64.  Furthermore,  a  dip  appears  at  the  approximate  wave  phase  velocity  in  the 
ion  distribution  curve  in  Fig.  12b  The  ion  phase  space  and  the  distribution  curve  are  presented 
in  Figs  I  3  a  and  13b  lor  17 /r,  =0.85.  The  corresponding  ion  phase  space  picture  For  /.  =4tt. 
and  NCi  =  128  is  given  in  Fig.  14,  showing  that  trapping  still  occurred  For  a  liner  mode  spacing. 
\i/  when  AA  was  reduced  to  hall',  i.e..  AA\/,  =0.07.  Figs  12.  13  and  14  show  the  dominant 
mode  at  saturation  to  he  AA,,  =  1/V2  ,  which  is  the  Fastest  growing  mode. 

In  general,  trapping  tends  to  Flatten  the  local  velocity  distribution,  centered  on  the  wave 
phase  velocity,  over  a  small  resonant  range  (  c,.,.—  +  r,,.J,.  F  From  Kq.  (8).  the  phase 

velocity  of  the  lower-hybrid  drill  instability. 


v,„  =  -  -  -f-ri)  .  (24) 

A  A-+A,;, 

are  well  separated  For  difFerent  mode  numbers  k,  and  are  much  less  than  the  ion  thermal  velo¬ 
city  For  i/  <  <  v:  Therefore,  each  mode  traps  independently  and  flattens  a  dilFerent  Fraction 
ol  the  ion  velocity  distribution  For  the  resonant  velocity  range  ()<r<r,/  Hence.  Hq.  120)  is 
replaced  by 


— 1 U= - \ — 7-  —  (25) 

4sv2n-  1+A-/A.;  [  r. 

For  the  saturation  level  of  each  discrete  mode  due  to  ion  trapping  Comparison  oF  the  simula¬ 
tion  saturation  Field  energy  spectrum  with  Fq,  (25)  is  given  in  Figs.  15a  and  15b  For  single¬ 
mode  and  many -mode  runs,  respectively.  The  co-existence  oF  the  multibeam  modes  causes  the 
large  error  bars  shown  in  Fig  15.  However  there  still  is  Fairly  good  agreement  between  simula¬ 
tion  and  theory  .  Both  simulation  data  and  the  theoretical  curves  For  small  A/A„  are  not  plotted 
because.  Following  F.q.  (X).  the  real  part  of  mode  Frequency  is  less  than  the  ion  cyclotron  Fre¬ 
quency  For  small  A/A,  .  and  the  lower-hybrid  drift  mode  turns  into  an  ion  drill  cyclotron 
mode,  our  theory  and  simulations  do  not  extend  into  that  region. 

I  rapping  was  also  observed  by  NVinske  and  l.iewer"  in  their  2d  particle  simulations  with  17 
large  than  r  For  >>\.  .  they  estimated  that  ion  trapping  required 
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(26) 
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anomalous  transport. 

f  or  time  i  =  n \i  >  mhh .  we  use  a  predictor-corrector  scheme  We  define  predicted  quan¬ 
tities  r(i  r  and  /,,,  as 


v.;,  =  v.l  1  and  =  T"  1  .  (29) 

B>  solving  Eq.  (15),  we  obtain  the  predicted  electric  held  and  the  held  energy.  We  then 

advance  the  particle  velocities  from  i"  1  :  to  »•// * 1 ' -  by  using  l'.',;  and  calculate  the  predicted  total 

ion  kinetic  energy  II ",  and  ion  momentum  , 

“ -  mi. x r"  '  1 :  (30) 

y  I 

P"r  =  (  p  1 :  +  A’;:'1  -  )/2  (3D 

\ 

where  P"'1  :  =  w  £»■",’ 1  :  .  Now  we  dehne  the  corrector  portion  of  the  algorithm.  Eq.  (27) 

i 

yields 


=  ‘V/!  +  <  P',1,- P"  )/niu, 

We  substitute  rj;  into  Eq.  (28)  and  obtain 


(32) 
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(33) 


We  then  substitute  r and  /,’  into  Eq.  (15)  and  solve  for  the  electric  held  P  and  the  held 


energy  « 


Using  this  method,  the  simulations  agree  with  the  linear  theory  at  small  wave  amplitudes 
for  various  values  of  nth.  Fig.  17  shows  a  held  energy  history  plot  using  time  varying  r,,  and 
/  and  / » y.  =0.85  and  ////; Ar=  1 00.  (The  results  are  not  sensitive  to  making  nth\t  smaller.) 
As  compared  to  Fig.  5b.  the  instability  now  saturates  at  a  much  lower  level  because  of  the 
current  relaxation  that  occurs.  By  assuming  that  the  system  is  stabilized  via  current  relaxation, 
the  saturation  level  of  the  fluctuation  spectrum  was  estimated  by  Davidson  el  al  h  at 


€,  =  — • 


|  inn,  vjj  2 


(34) 


2  2(1  +ca,',,  I<d~,  ) 

(he  corresponding  changes  in  the  ion  temperature  and  electron  temperature  after  saturation 


were  given  by 
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lor  /  /  /  <  <  I  and  » /  / 1 ...  «  I .  For  our  simulation  parameters,  A  7 ,11  and  A///  are  about 
10  '  .  and  are  too  small  to  be  observed.  Comparison  of  simulation  results  lor  time  varying  i ,, 
and  /  with  that  for  constant  i„,  and  }\  is  given  in  Fig  10.  showing  that  current  relaxation  sta¬ 
bilizes  the  instability  when  i,(i  and  7,  are  dynamic,  otherwise,  ion  trapping  stabili/es  the  insta¬ 
bility  For  cases  in  which  r,  <<i  and  i,,  is  allowed  to  relax  in  lime,  the  saturation  level  of 
t./nl  is  so  small  that  it  is  masked  by  the  thermal  fluctuations,  whose  levels  are  equal  to 
~  10  '  in  our  simulations.  Thus,  no  data  appears  on  the  curve  labeled  lb)  in  Fig.  10  for 
\,  i  <0.8  .  Note  that  m  a  finite  beta  plasma,  the  particle  drifts  and  magnetic  field  are  coupled 
and.  consequently,  the  anomalous  dissipation  of  the  particle  drift  energy  and  magnetic  energy 
are  also  linked.  Therefore,  the  magnetic  held  effectively  acts  as  a  free  energy  source  to  drive 
ihe  instability,  and  current  relaxation  will  not  occur1*  . 


\  .  CONCT.l  SIONS 

One-dimensional  particle-hybrid  simulations  of  the  lower-hybrid  drift  instability  in  the  low 
drift  velocity  regime  have  been  described  in  this  paper.  At  small  amplitudes,  the  simulations 
agree  very  well  with  the  linear  theory.  Saturation  mechanisms  have  been  investigated  for  cross- 
held  drill  velocities  less  than  the  ion  thermal  velocity  It  has  been  shown  that  the  lower-hybrid 
drift  instability  is  stabilized  by  ion  trapping  if  c.oss-field  drift  velocities  are  kept  constant  in 
time  because  of  external  circuits,  for  example.  An  analytic  formula,  which  agrees  very  well 
with  simulation  results,  lor  the  saturated  field  energy  spectrum  m  k  space  for  ion  trapping  has 
been  given  in  See  III  A  If  cross-held  drift  velocities  arc  allowed  to  vary  in  time,  growth  ol  the 


instability  will  be  stabilized  through  current  relaxation 
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to,,.  /- 600  lor  \f  =  kL/2n—5.  parameters  of  F  ig  2  The  sharp  peak  ai  lu/tu,.  =0.1.16 


corresponds  in  the  lower-hvhrid  drill  instability 
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Fleet  ric  field  energy  history  plots  for  parameters  of  Fig.  2,  except  for  (a)  t,  /i  =0.57  . 
the  high  frequencies  are  due  to  the  multibeam  instability  and  low  frequency  is  due  to 
the  LI  11)1.  and  <b)  y,/v  =0.85  ;  here  the  multibeaming  effects  are  much  smaller  at 
saturation,  which  is  much  higher  than  in  (a) 


ower-hvbrid  drill 


f  ig  8  Lower-h>hrid  drill  dispersion  curves  for  17 /r,  =1.5556.  (;i> 


because  that  only  one  of  every  5  panicles  was  pinned  The  space  between  beams  i 


f  ig.  10  Comparison  of  saturation  field  energies  as  functions  of  \-,/r, ..  (a)  For  constant  i,,( 
and  /, ,  the  lower-hybrid  drift  instability  is  saturated  by  ion  trapping,  (b)  For  time 
varying  \llt  and  /, .  saturation  of  the  instability  is  achieved  via  current  relaxation. 
The  simulation  data  (dots  in  (a)  and  crosses  in  (b)  )  are  in  better  agreement  with  the 
theory  in  (a)  than  in  (hi 


F  ig.  1 1  Comparison  of  simulation  trapping  fre 
theory.  Hq.  (23).  oii=2tt/T,  as  shown 


city  distribution  function  after  saturation  ls  the  phase  velocity  of  the  most 
unstable  mode 


vc loci t\  distribution  (unction  after  saturation  r,„  is  the  phase  \elociH  of  the  most 


